First Test of Lorentz Violation with a Reactor-based Antineutrino Experiment 
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We present a search for Lorentz violation with 8249 candidate electron antineutrino events taken 
by the Double Chooz experiment in 227.9 live days of running. This analysis, featuring a search for 
a sidereal time dependence of the events, is the first test of Lorentz invariance using a reactor-based 
antineutrino source. No sidereal variation is present in the data and the disappearance results are 
consistent with sidereal time independent oscillations. Under the Standard-Model Extension (SME), 
we set the first limits on fourteen Lorentz violating coefficients associated with transitions between 
electron and tau flavor, and set two competitive limits associated with transitions between electron 
and muon flavor. 



Recently, we reported evidence of electron antineu- 
trino disappearance with the Double Chooz far detec- 
tor, 1050 m away from two 4.25 GW reactor cores [l|, Q, 
which generally is interpreted in terms of mass-induced 
neutrino oscillations. A path to more exotic physics be- 
yond the Standard Model (SM) may be gained by care- 
fully examining the oscillation behavior. In particular, 
the collected electron antineutrino sample also provides 
an opportunity to search for the violation of Lorentz in- 
variance. 

Lorentz invariance requires that the behavior of a par- 
ticle is independent of its direction or boost velocity. The 
as-yet-unseen violation of this principle is predicted to 
occur at the Planck scale and is especially interesting as 
it can occur dynamically via spontaneous Lorentz sym- 
metry breaking Q. The process of neutrino oscillation, 
in which a neutrino of one flavor transforms into another 
flavor after traveling a distance, is due to interference 
between the slightly different Hamiltonian eigenstates of 
the propagating particle. The experimental observable, 
oscillation probability, is therefore quite sensitive to small 
couplings between neutrinos and a possible Lorentz vio- 
lating field. 

Testing Lorentz violation with the natural interferom- 
eter of neutrino oscillation has been done in several ex- 
periments, including MINOS H, IceCube @, LSND @, 
and MiniBooNE 0|. These tests all fall under the for- 
malism of a search within the Standard- Model Extension 
(SME) [8.]. In this letter, we describe the first search for 
Lorentz violation using reactor antineutrinos. 

In the SME, all possible types of Lorentz violation are 
added to the SM Lagrangian. Here, we limit ourselves 
to the renormalizable sector (referred to as the minimal 
SME). For Lorentz violating neutrino oscillation, the ef- 
fective Hamiltonian is written as Q 

(K s u ~ + |[(oi)^ - (c L r P , P „} ab , (i) 

where E and p^ are the energy and 4-momentum of the 
neutrino and (m 2 ) ab refers to the neutrino mass in the 
flavor basis represented by a and b. The CPT-odd coeffi- 
cient (at)^ switches sign for antineutrinos and violates 



both Lorentz and CPT symmetry, while the CPT-even 
coefficient (cl)^ violates Lorentz but maintains CPT. 
Both are vector and tensor and consist of direction in- 
dependent parts [{a L )l b , (a L )f b , (c L )^ T , (ci)^f, and 
(c L )ab Z ] and direction dependent parts [(a L )* b , {a L ) Y b , 
(cL)l b X , {c L ) T ab Y , {cL)* b X , {c L )^ b Y , (c L )X b z , {c L ) Y ab Y , and 
{ c Llab\ m * ne Sun-centered coordinate system (repre- 
sented by the superscripts). A measured non-zero direc- 
tion dependent component would be clear evidence of an 
anisotropy in the universe and Lorentz violation. As dis- 
cussed later, no evidence for Lorentz violation has been 
found and our goal is therefore to set limits on these co- 
efficients. We note that the known neutrino mass term in 
the flavor basis in this formalism is neglected in order to 
follow a conservative approach when setting these limits. 

Sidereal time is based on the Earth's orientation rel- 
ative to the fixed stars. The unambiguous signature of 
Lorentz violation is a sidereal modulation of an experi- 
mental observable such as neutrino oscillation probabil- 
ity. A sidereal variation is expected for an experiment 
moving in a fixed Lorentz violating field with the rota- 
tion of the Earth. We probe this field by searching for 
such a dependence among the collected electron antineu- 
trino events. The antineutrino vector is set using the an- 
tineutrino source and the location of the detector. The 
location of the source is taken to be a weighted point in 
between the two cores, 6° apart relative to the detector, 
representative of the number of antineutrinos expected 
from each during the physics run. 

The data set used for this analysis was obtained with 
the Double Chooz experiment between April 13, 2011 
and March 15, 2012. Electron antineutrinos interact in 
the detector via the inverse beta decay (IBD) process, 
D e + p — >• e + + n. IBD events produce a distinct dou- 
ble coincidence signature from the prompt positron sig- 
nal followed by neutron capture 30 /is (mean time) later. 
The "inner detector," composed of three concentric cylin- 
drical regions separated by acrylic, is used to observe and 
efficiently reconstruct these events as well as mitigate 
background. The innermost 10 m 3 cylinder contains 1 g/1 
gadolinium-doped scintillator and forms the antineutrino 
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target. Surrounding this is the "gamma catcher," de- 
signed to detect gamma-rays escaping the target volume. 
The gamma catcher volume is then enveloped by a non- 
scintillating oil buffer in which 390 10 inch PMTs are 
immersed. The inner detector is surrounded by a steel 
vessel that forms an optically isolated outer cylinder filled 
with scintillator. This "inner veto" , along with an "outer 
veto" mounted above it and 15 cm of shielding steel, is 
used to reject cosmic ray events. 

The antineutrino sample and event selection criteria 
are identical to those used for the disappearance analy- 
sis reported in Ref. 0] . The data consist of 8249 IBD 
candidates, collected with 227.9 live days and 33.7 GW- 
ton- years exposure. There are 497 background events 
expected in this sample. The background is mainly com- 
posed of (1) cosmogenic radioisotopes, such as 8 He and 
9 Li, which decay via the emission of (3n, (2) cosmogenic 
stopping muons as well as fast neutrons that interact 
multiple times in the inner detector, and (3) acciden- 
tal coincidence of a radioactivity-induced prompt signal 
followed by a neutron-like event. Background event rate 
as a function of sidereal time is treated as a constant. As 
the dominant background contributions to the Double 
Chooz analyses arise as the result of cosmic ray muons, 
we study the time dependence of muon veto rate in order 
to justify this assumption. The maximum variation in 
muon veto rate as a function of sidereal time is about 
0.5%. A background variation in time at this level would 
create a maximum variation in disappearance probability 
of -0.03%. 

The background-subtracted IBD candidate sample is 
directly compared to the Monte Carlo (MC) expectation 
in order to probe a possible sidereal time dependence. 
The unoscillated MC expectation is based on the IBD 
cross section, the reactor flux prediction, the detector 
response, and the number of protons in the detection 
volume. The expectation is formed from each of these 
variables on a run-by-run basis, with each physics run 
lasting approximately one hour. We note that the ther- 
mal power of each core is estimated in <1 minute time 
intervals and the uncertainty on the total power is 0.5%. 
The reactor flux prediction uses extensive input from the 
Chooz reactor facility and Electricite de France (EPF). 
The quality of the code has been benchmarked [lOj and 
compared to EDF assembly simulations. The P e spec- 
trum is taken directly from Refs. Ill 1 1 211 and is normal- 
ized to the Bugey4 rate measurement [13j . The analysis 
input information, shown in Fig. [TJ is assigned to one of 
24 bins between and 23.934 hours (one sidereal day). 
A MC expectation event weight is split up between the 
relevant time bins based on the time and length of the 
run while a data event is placed in a bin based on its 
DAQ timcstamp. 

A number of sources of systematic uncertainty are con- 
sidered. These include those associated with the back- 
ground prediction, the detector and detector response, 
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FIG. 1. The background subtracted data and MC expecta- 
tion IBD event rates as a function of sidereal time. The MC 
expectation assumes no antineutrino disappearance. Total er- 
rors (statistical and systematic) are shown on the data points. 



and the reactor flux (normalization and shape). The re- 
actor flux and detector operations are both weak func- 
tions of solar time due to human activity as the cores 
turned on/off multiple times during the run and detec- 
tor calibrations are generally done during the daylight 
hours. Day-night effects are well accounted for in the MC 
prediction and correlated uncertainties associated with 
these are included in a covariance matrix, fully describ- 
ing the predicted statistical and systematic errors. The 
3.93 minute/day difference between sidereal and solar 
time, compounded over the ~1 year physics run, largely 
removes any potential for an unaccounted modulation 
in sidereal time associated with small modulations re- 
lated to solar dependence. The detector and background 
prediction uncertainties arc considered uncorrelatcd with 
each other and fully correlated in sidereal time. A thor- 
ough explanation of the various uncertainties and their 
determination can be found in Ref. 0], noting that cor- 
relations in time (as opposed to antineutrino energy in 
the reference) are most important here. The total frac- 
tional uncertainty with respect to the MC expectation is 
2.9%. The statistical uncertainty contributes at the level 
of 1.1% and systematic uncertainties are led by the reac- 
tor flux and detector response (1.7%) and the background 
prediction (1.7%). 

In the three active flavor neutrino oscillation frame- 
work, the v e — > v e probability can be written as a 
function of D e — > and v e — > v T (Pp e ^p e = 1 — 
P^^-Pp^pJ. Under the SME, both P^^ and 
Pp^v are written as functions of five free parame- 
ters 
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P*„ 



|(^off)e/i| 2 i 2 \(h e s)ef\ 2 L 2 



(He) 



{hey 



L 2 



- - l I (C)efi + (A)ep sina; e r ffi + (A c ) sp , cosw e T e 

+ (B s )ep,sm2u;®T @ + (B c ) Sfi cos2w e T e | 2 

+ | (C) gf + (A s ) S f sino; ffi T ffi + (A c )sf cosw ffi T e 

+ {B a ) S r sin 2w e Te + (8 c )er cos 2w e T e I 2 ] ( 2 ) 
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The disappearance probability is a function of sidereal 
time T®, sidereal frequency w® [2n/ 86 164.1 base- 
line L, and ten amplitudes (parameters). The param- 
eters themselves are composed of the Lorentz violating 
coefficients introduced in Eq. [TJ antineutrino energy, and 
the antineutrino-source-to-detector vector. We aim to re- 
duce this equation since there are too many parameters 
for a realistic fit and measurement extraction. Ideally, 
this reduction proceeds without any assumptions, in a 
model independent way. 

Double Chooz's maximum sensitivity to the CPT- 
odd and CPT-even SME coefficients is on the order of 
~10 -20 GeV and ~10~ 18 , respectively, determined by 
considering the maximum oscillation condition in the ef- 
fective Hamiltonian (Eq. H). The MINOS near detec- 
tor Q and MiniBooNE [tJ measurements both place sig- 
nificantly better limits on all CPT-even coefficients, at 
the level of ^-T0~ 21 and ^10 -20 , respectively. The ten 
relevant SME coefficients are therefore set to zero, corre- 
sponding to the removal of two parameters, (B s )gp. and 
{13c)efn from Eq. [2l It is now difficult to remove more 
parameters in a model independent way and we cannot 
reduce Eq. [2] further using existing measurements. We 
therefore study two different sets of assumptions. 

The assumption that all Lorentz violating oscillations 
occur in electron antineutrino to tau antineutrino transi- 
tions (P Pc ^p ti = 0, P Pe ^p T ^ 0) is studied with the "e — r 
fit": 

L 2 

Pp e ->P e ^ 1 - 7^32 [ I ( C )sf + (As)er sillCJ®T e 

+ {A c )ef cosw e r e + (B s )gr sin2w ffi T ffi 
+ (23 c ) gf cos2w e Te| 2 ] (3) 

The five free parameters [(C) S f , (A s ) S f, (A)g?, (B 8 )et, 
and (B c )ef] themselves contain fourteen of the e—r sector 
SME coefficients introduced in Eq. [T] 

The second model is based on the assumption that 
all Lorentz violating oscillations occur in electron an- 
tineutrino to muon antineutrino transitions (P Pe _> 5 7^ 0, 
Pv a -n? T = 0) and is referred to as the "e — /x fit": 

Pp c ^p c ^ 1 - [ I ( C W* + (A)ejs sinw e T ffi 

+(A) g pCosw e r e | 2 ] (4) 

This equation has only three free parameters [(C) gp,, 
(A s )gji, and (Ac)ep\, as the MINOS and MiniBooNE 
constraints have removed the CPT-even coefficients, and 
contains four e — 11 sector SME coefficients. The C pa- 
rameter in Eqs. [3] and HI contains sidereal time indepen- 
dent SME coefficients. This term can affect both shape 
and normalization. We note that each of these two mod- 
els considers disappearance in only one channel while 
the complete formula (Eq. [2]) contains contributions from 
both. However, the limits reported tend to be more con- 
servative than if both channels were considered simulta- 
neously. 



The SME parameters in Eqs. [3] and [4] are extracted 
using the MC expectation and background subtracted 
data (Fig. [T]), the total error matrix including statisti- 
cal and correlated systematic contributions, and a least 
squares fitting technique. The least squares estimator is 
minimized in order to find the best fit (BF) among the 
parameter combinations. 

The e — r and e — /x BF sidereal time results are shown 
in Fig. [2] The BF results for both fits are dominated 
by the sidereal time independent terms, (C)st and (C)sp- 
We examine the significance of the results below. 
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FIG. 2. The electron antineutrino disappearance probability 
as a function of sidereal time, overlaid with the best fit e — fi 
(x 2 /ndf = 28.8/21) and e - r ( X 2 /ndf = 27.7/19) curves. A 
one parameter flat fit to the distribution yields a best fit with 
X 2 /ndf = 30.6/23. 

The flatness of the sidereal time distribution is ana- 
lyzed using a frequentist approach. A large sample of 
randomized pseudoexperiments, based on the MC expec- 
tation and the total error matrix and with an injected 
sidereal time independent ( "flat" ) disappearance, is gen- 
erated in order to determine the fraction of samples that 
present a more or less flat solution than the one found 
here. We introduce a normalization factor of 91.8%, con- 
sistent with the counting-only disappearance probability, 
in order to ensure that we are testing the null hypothesis 
that there is no sidereal time dependence rather than the 
null hypothesis that there is no antineutrino disappear- 
ance. The Ax 2 is defined as the minimum \ 2 from the 
flat hypothesis minus the minimum \ 2 from each e — r 
or e — n fit. This frequentist study shows that 60.0% 
(41.8%) of pseudoexperiments have a larger A\ 2 than 
the real data and that the e — r (e — n) results are consis- 
tent with sidereal time independent oscillations. In the 
absence of a sidereal dependence, we proceed to set limits 
on the relevant time dependent SME coefficients. 

Limits on the SME coefficients and uncertainties asso- 
ciated with the BF parameters are determined by con- 
structing a five and three dimensional parameter space, 
corresponding to the e — r and e — ll fits, respectively. By 
assuming the minimum of the least squares fit estimator 
follows a x 2 distribution, a \a (2a) hyper-volume can be 
defined as the region enclosed by the constant \ 2 hyper- 
surface with minimum \ 2 plus 5.9 (11.3) for the e — r fit, 



5 





BF parameter 


1(7 limit 


SME coefficients combination 




(10- 20 GeV) 






(C) g f 


5.8 


7.8 
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(C)ep 
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[-(a i )^-0.29(a i )f,J 
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-0.4 ±0.7 


1.9 


[-0.91(a i )f M ±0.29(a i )r M ] 


(Ac)ejl 


0.5 ±0.8 
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TABLE I. A summary of the e — r and e — [i Lorentz violation measurements in terms of the best fit (BF) parameters and the 
corresponding combinations of Standard-Model Extension coefficients. The average antineutrino energy "E" is 4.2 x 10 -3 GeV. 
Note that the e — r and e — ft fit functions are not linear with the fit parameters and la errors do not scale to the 2a ones. 



and 3.5 (8.0) for the e — fx fit. These criteria are tested 
by using a sample of pseudoexperiments with an injected 
signal based on the BF. That is, each pseudoexperiment 
sample is convolved with the BF oscillation probability 
equation. A new fit is then performed and the BF param- 
eters are tallied. We find that the above choices for la 
and 2a hyper-surfaces enclose 70.1% (68.8%) and 94.5% 
(94.7%) of BF points for the e — r (e — fi) fits and that our 
allowed regions are valid. Note that we have considered 
only half of the parameter space in this procedure and 
that the sign reversed BF parameters are equally valid. 

The results are summarized in Table HI The BF val- 
ues from both the e — r and e — fj, fits are shown along 
with la errors and 2a upper limits, when applicable. The 
BF errors are generally asymmetric; however, the larger 
of the two-sided error bar is reported. Correlations be- 
tween parameters and multiple connected solutions make 
it impossible to extract meaningful la errors on the BF 
e — r parameters. The combination of SME coefficients 
associated with each measured parameter is also shown 
in the table. All e — r parameter limits as well as the 
sidereal time dependent e — /i limits are on the order of 
~10~ 20 GeV for CPT-odd coefficients and ~10~ 17 for 
CPT-even coefficients. 

Although every measured sidereal time dependent pa- 
rameter is consistent with zero, the time independent 
parameter (C) g p is non-zero at the 2.1a level. We note 
that a normalization-only fit ~ 1 — ((C)g f + 

(C)| p )] yields (C% + (C)% = (34.2 ±9.2) x lO" 40 GeV 2 . 
This disappearance is consistent with the rate-only #13 
measurement in Ref. Q. With current precision, time 
independent Lorentz violating effects cannot be distin- 
guished from mass and #13 induced oscillations. Sepa- 
rating the two effects may be possible with future high 
statistics data and spectral information, however. The 
disappearance observed can generally be interpreted as 
due to neutrino mass and #13 in the three flavor neutrino 
oscillation framework. 

There are a number of alternative neutrino oscillation 
models motivated by Lorentz violation 15- 1^. These 



models neglect sidereal modulations by assuming that 
any such variations are averaged out or that the probabil- 
ity of oscillation is governed by time independent terms 
only. The models focus on reproducing the global ob- 
served energy and baseline dependence of neutrino os- 
cillations. Interestingly, however, none of the models 
predict the observed antineutrino oscillations at Double 
Chooz's energy ((E}=4.2 McV) and baseline (1050 m). 
That is, the measured disappearance conflicts with these 
models. This may be an additional reason to interpret 
the time independent disappearance observed as due to 
neutrino mass and non-zero #13, rather than time inde- 
pendent Lorentz violation. 

We have analyzed the sidereal time dependence of Dou- 
ble Chooz's electron antineutrino candidates as a probe 
of Lorentz violation. With no observed modulation, we 
set the first limits on fourteen of the SME coefficients in 
the e — r sector, and set competitive limits on two e — fi 
sector coefficients. Competitive limits may also be pro- 
vided by other reactor antineutrino experiments in the 
future 2(], 21 1, although Double Chooz features a com- 
paratively simple antincutrino-source-to-detector vector. 
With the addition of this work amongst the world's data, 
sidereal variation tests with neutrino oscillation exper- 
iments have been performed with all active oscillation 
channels. In the future, astrophysical neutrinos |22| may 
improve sensitivity to Lorentz violation by many orders 
of magnitude compared to what is possible for terrestrial- 
based neutrino experiments. 
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